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I. INTRODUCTION 



Soliton solutions to nonlinear partial differential equations (PDEs) in (1+1) dimensions 
have been of great interest ever since the very discovery of completely integrable PDEs. 
There is a wide range of literature concerning integrable nonlinear PDEs and their soliton 
solutions (see, for instance, Refs. |l|-§l and the references therein). Much is known about the 
behavior of solitons and their interactions in various integrable systems (soliton scattering, 
breather solutions, soliton bound states, etc). Such knowledge is very valuable not only 
for the underlying integrable systems, but also for nearly integrable systems which can be 
studied analytically by soliton perturbation theories. 

It is an astonishing fact, however, that notwithstanding the almost three decades of 
advances in the study of soliton dynamics, still there are substantial gaps in our knowledge 
of soliton solutions in integrable nonlinear PDEs. Indeed, the reader familiar with the 
basics of the inverse scattering transform method knows that it is poles of the reflection 
coefficient (or, in modern terms, zeros of the Riemann-Hilbert problem) that give rise to 
the soliton solutions. The soliton solutions are usually derived by using one of the several 
well-known techniques, such as the dressing method or the Riemann-Hilbert problem 

approach However, in most publications (except, to our knowledge, Refs. |f^-[13[|, which 
are discussed below) only soliton solutions from simple poles are considered. It is usually 
assumed that a multiple-pole solution can be obtained in a straightforward way by coalescing 



several distinct poles (see, for instance, Ref. PJT^) which describe multisoliton solutions. 
This would be indeed the case if such coalescing were a regular limit. However, this limit is 
obviously a singular one. Indeed, the soliton dressing matrix corresponding to a multisoliton 
solution is a rational matrix function which has distinct simple poles, while the coalescing 
procedure must produce multiple poles. Obviously, a more careful examination of this issue 
is necessary. This is the main subject of the present paper. 

Soliton solutions corresponding to multiple poles, i.e., the higher-order solitons, have 
been investigated in the literature before. A soliton solution to the nonlinear Schrodinger 
(NLS) equation corresponding to a double pole was first given in Ref. [0] but without much 
analysis. The double- and triple-pole soliton solutions to the KdV equation were exam- 
ined in Ref. and the general A^-pole soliton solution to the sine-Gordon equation was 
extensively studied in Ref. using the associated Gelfand-Levitan-Marchenko equation. In 
Refs. ||T^,[n]], higher-order soliton solutions to the NLS equation were studied by employing 
the dressing method. Finally, in |jl2|,[13|, higher order solitons in the Kadomtsev-Petviashvili I 
equation were derived by the inverse scattering method. 

In this article, we consider the higher-order zeros of the Riemann-Hilbert problem for the 
A^-wave system and study the corresponding soliton matrices. We derive the soliton dressing 
matrices for the simplest class of higher-order zeros - the elementary higher-order zeros. We 
call the n-th order zero ki of the Riemann-Hilbert problem elementary if the soliton matrix 
evaluated at ki has only one vector in the kernel, i.e. the geometric multiplicity of the zero 
is 1 (see Definition 1 in Sec. 3 for more details). The corresponding higher-order soliton 
solutions to the iV-wave system are determined. Then we apply our theory to the physically 
important three- wave interaction model and derive the simplest higher-order soliton solution. 
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In the generic case, this solution describes the breakup of a higher-order pumping wave into 
two higher-order elementary waves, and the reverse process. In non-generic cases, this 
solution could describe {i) the merger of a pumping sech wave and an elementary sech wave 
into two elementary waves (one sech and the other one higher-order); {ii) the breakup of a 
higher-order pumping wave into three sech waves — one pumping wave and two elementary 
waves. The higher-order soliton solution could also reproduce fundamental soliton solutions 
special case. 

In general, one needs to consider zeros with the geometric multiplicities taking values 
from 1 to — 1, where is the matrix dimension of the Riemann-Hilbert problem. The 
present work is the first step towards the solution of this general case. The point is that the 
soliton matrices derived here for the elementary zeros provide the building blocks for the 
most general case. We plan to address the general problem in the next paper. Thus, for the 
3-wave interaction model, the higher-order soliton solutions we derived in this paper (which 
correspond to elementary higher-order zeros) may not be the most general higher-order 
soliton solutions in this wave system. 

It is noted that the soliton dressing matrices for the higher-order soliton solutions were 
already a subject of interest in Refs. [|TO|,|lT|, where an ansatz for the soliton matrices was 
proposed for the 2x2 Zakharov-Shabat spectral problem and expressions for the higher-order 
soliton solutions of the nonlinear Schrodinger equation were obtained. In the present paper, 
we study higher-order solitons in the A^-wave system. Even though the evolution equations 
considered in |lO|,|ll| and this article are different, their Riemann-Hilbert formulations have 
a lot in common. In this paper, we place the emphasis on the derivation of the soliton 



matrices for the higher-order zeros and fill some gaps in the approach of Ref. |Tl]]. We 
show that the ansatz of |Tl| is precisely the soliton matrix for an elementary higher-order 
zero in a X spectral problem. For the 2x2 Zakharov-Shabat spectral problem, any 
higher-order zero is elementary. Thus the ansatz of |TI| is the general soliton matrix of the 



Zakharov-Shabat problem. Consequently, higher-order soliton solutions obtained in |10,11 



are the general higher-order soliton solutions in the nonlinear Schrodinger equation. But in 
a N X N spectral problem with A^ > 2, a higher-order zero is non-elementary in general. In 
that case, the ansatz of |TT[ will not be the general form of soliton matrices. In this paper, 
we also uncover the invariance property of the soliton matrix for a higher-order zero, which 
is necessary for the vector-parametrization of the soliton matrix to be self-consistent. The 
invariance property obtained in [|11[] is shown to be just a special case. Lastly, we point 
out that our derivation of higher-order solitons in the A^-wave system is made for general 
dispersion laws and arbitrary matrix dimensions. In addition, our results can be generalized 
to the most general case of non-elementary zeros of the N x N spectral problem. 

The paper is organized as follows. A summary on the Riemann-Hilbert problem is placed 
in section ||. Section |T| is the central section of the paper. There we present the theory of 
soliton matrices corresponding to the higher-order zeros of the Riemann-Hilbert problem. 
For an elementary higher-order zero of the N x N Riemann-Hilbert problem (see definition 
in the text), the general soliton matrix is derived. This soliton matrix is similar to the ansatz 
as proposed in |TT|] . In section |V|, the theory is applied to derive the simplest higher-order 
soliton solution in the three-wave interaction model. 
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II. THE RIEMANN-HILBERT PROBLEM APPROACH: SUMMARY 



The integrable nonlinear PDEs in 1+1 dimensions are associated with the matrix 
Riemann-Hilbert problem (consult, for instance, Refs. |]T]-^,|T^-|2^). The matrix Riemann- 
Hilbert problem (below we work in the space of N x N matrices) is the problem of finding 
the holomorphic factorization, denoted below by ^+{k) and ^Z^{k), in the complex plane 
of a nondegenerate matrix function G{k) given on an oriented curve 7: 

^Z^ik,x,t)^+{k,x,t) = G{k,x,t) 

= exp[-A{k)x-n{k)t]G{k,0,0)exp[A{k)x + n{k)t], ke-f. (2.1) 

Here the matrix functions $+(/c) and $I^(A;) are holomorphic in the two complementary 
domains of the complex fc-plane: (7+ to the left and C_ to the right from the curve 7, 
respectively. The matrices A(fc) and Q{k) are called the dispersion laws. In this paper, we 
require the dispersion laws to be diagonal (we have accounted for this on the right-hand 
side of equation ( p.l|) by writing the explicit (x, t)-dependence in the form of an exponent). 
The Riemann-Hilbert problem requires an appropriate normalization condition. Usually 
the curve 7 contains the infinite point /c = 00 of the complex plane and the normalization 
condition is formulated as 

$±(/c,x, t) ^ /, as k ^ 00. (2.2) 

This normalization condition is called the canonical normalization. Setting the normalization 
condition to an arbitrary nondegenerate matrix function S{x, t) leads to the gauge equivalent 
integrable nonlinear PDE, e.g., the Landau-Lifshitz equation in the case of the NLS equation 
1^. Obviously, the new solution ^±{k,x,t) to the Riemann-Hilbert problem, normalized to 
S{x,t), is related to the canonical solution by the following transformation 

$±(fc, X, t) = S{x, t)^ik, X, t). (2.3) 

Thus, without any loss of generality, we confine ourselves to the Riemann-Hilbert problem 
under the canonical normalization. 

For physically applicable nonlinear PDEs the Riemann-Hilbert problem possesses the 
involution properties, which reduce the number of the dependent variables (complex fields). 
The A^-wave interaction model admits the following involution property of the associated 
Riemann-Hilbert problem 

^l{k) = ^Z\k), k = k*. (2.4) 

Here the superscript "f" represents the Hermitian conjugate, and "*" the complex conjugate. 
However, our approach can be trivially extended to the general case without such involution. 
To keep our treatment general, we will use the overlined quantities where applicable. The 
reduction to the involution is then done by associating the overline with the Hermitian 
conjugation in the case of vectors and matrices and with the complex conjugation in the 
case of scalar quantities. 
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To solve the Cauchy problem for the integrable nonlinear PDE posed on the whole axis 
X, one usually constructs the associated Riemann-Hilbert problem starting with the linear 
spectral equation 

(9^$(fc, X, t) = (^{k, X, t)A{k) + U{k, X, t)$(fc, X, t), (2.5) 

whereas the t-dependence is given by a similar equation 

dt^{k, X, t) = ^{k, X, t)n{k) + V{k, x, t)^{k, x, t). (2.6) 

The nonlinear integrable PDE corresponds to the compatibility condition of the system ( p.5| ) 
and (12:^ ): 



dtU - d,V + [U, V] = 0. 



(2.7) 



The essence of the approach based on the Riemann-Hilbert problem lies in the fact that 
the evolution governed by the complicated nonlinear PDE ( p.7|) is mapped to the evolution 
of the spectral data given by simpler equations such as ( p.l|) and ( p.21| ). For details, consult 
Refs. iilMl. 

Let the evolution equations for the spectral data be given. In our case, these are equation 
( ^yil) for G and equation (p.21|) (see below) for the discrete data. Then the matrices U{k, x, t) 
and V{k,x,t) describing the evolution of $± can be retrieved from the Riemann-Hilbert 
problem. In our case, the potentials U{k, x, t) and V{k, x, t) are completely determined 
by the (diagonal) dispersion laws A(A;) and VL{k) and the Riemann-Hilbert solution $ = 
$±(A;,x,t). Indeed, let us assume that the dispersion laws are polynomial functions, i.e.. 



Jl J2 

A(fc) = Ajk\ n{k) = ^ Bjk^. 

j=0 j=0 



Then using similar arguments as in Ref. [T^ we get: 

?7 = -p{$A$-^}, v = -v{m^-^}. 

Here the matrix function $(A;) is expanded into the asymptotic series, 

$(A;) = / + A;"^<I>W +A;"2<I>(2) + ..., ^ oo, 

and the operator V cuts out the polynomial asymptotics of its argument as k 
important property of matrices U and V is that 

TrU{k, X, t) = -TiA{k), TiV{k, x, t) = -Trfi(A;), 



(2.8) 



(2.9) 



oo. An 



(2.10) 



which evidently follows from equation ( |2.9| ). Below, let us consider the three- wave interaction 
system as an example [3,p5H27||. Set = 3, 



^ai ^ 



A(A;) = ikA, A 



a2 
ag 



n{k) = ikB, B 



^ hi ^ 

62 
63 



[2.11) 



where aj and bj are real with the elements of A being ordered: ai > 02 > 03. From equation 

we get 



U = -A{k)+i[A,<l>^% V^ = -1](A;) (2.12) 

Setting 



u 



1 = V^i - 0.2^12^ U2 = y/a2 - a3^23^ ^3 = V^i " ct^^is ^ (2-13) 



assuming the involution (|2.4| ), and using equation (|2.12| ) in (|2.7|) we get the three- wave 
system: 

dtUi + vidxUi + ieu2U3 = 0, (2.14a) 

c?iM2 + V2dxU2 + ieuiUz = 0, (2.14b) 

dtU3 + vadxUs + ieuiU2 = 0. (2.14c) 



Here 



b2-bi bs-b2 bs-bi 
vi = , V2 = , V3 = , (2.15) 



_ 0162 - a2&i + a2&3 - «3&2 + 03^1 - «1&3 
[(ai - a2)(a2 - a3)(ai - ^3)]^/^ 

The group velocities satisfy the following condition 

V2 - f 3 ai - 02 



(2.16) 



vi - f 3 02 - as 



< 0. (2.17) 



The three-wave system ( p.l4| ) can be interpreted physically. It describes the interaction 
of three wave packets with complex envelopes ui, U2 and M3 in a medium with quadratic 
nonlinearity. 

It is often desirable to relate the inverse-scattering parameters aj and bj {j = 1,2,3) to 
the physical parameters e and Vj (j = 1, 2, 3). This relation can be easily found from (p.l5|) 
and (|2l6| ) as 



ai-a2 = - -, 02-03 = 7 TT r- (2.18) 

(^1 - ^2)(^^l - ^^3) (^1 - t^2)(^^3 - ^^2) 

The other parameters ai — and bj {j = 1,2,3) can be readily obtained from equations 
( ^.18| ) and ( |2.15| ). Note that the inverse-scattering parameters are not uniquely determined. 
In fact, one of aj and one of bj {j = 1,2,3) are free parameters. It is an invariance in 
the inverse-scattering formulation of the 3-wave system and it does not affect the physical 
solution in any way. 

In general, the Riemann-Hilbert problem ( p.l| )-( pl2D has multiple solutions. Different 



solutions are related to each other by the rational matrix functions T{k) (which also depend 
on the variables x and t) P,|3|j5|,|^, |l^] : 
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(2.19) 



The rational matrix r(A;) must satisfy the canonical normalization condition: T{k) — > / for 
k oo and must have poles only in C_ (the inverse function r~^(A;) then has poles in C+ 
only). Such a rational matrix r(A;) will be called the soliton matrix below, since it gives the 
soliton part of the solution to the integrable nonlinear PDE. 

To specify a unique solution to the Riemann-Hilbert problem the set of the Riemann- 
Hilbert data must be given. These data are also called the spectral data. The full set of the 
spectral data comprises the matrix G{k, x, t) on the right-hand side of equation (|2.1| ) and the 
appropriate discrete data related to the zeros of det $+(fc) and det $I^(/c). We will confine 
ourselves to the case of the Riemann-Hilbert problem with zero index, i.e., when det $+(/c) 
and det$I^(A;) have equal number of zeros (counting the multiplicity). For instance, in 
the case of involution ( p.4|) the Riemann-Hilbert problem has zero index because the zeros 
appear in complex conjugate pairs: kj = k*. It is known p9|-pl| (see also Ref. [Q) that in 
the generic case the spectral data include simple (distinct) zeros fci, . . . , /c„ of det ^+{k) and 
ki, . . . ,kn of det $l^(/c), in their holomorphicity domains, and the null vectors . . . , 
and {vi\, . . . , {Vn\ from the respective kernels: 

<^+{kj)\v,) = 0, {vj\^Z\kj) = 0. (2.20) 

Using the property ( p.lO|) one can verify that the zeros do not depend on the variables 
X and t. The (x, t)-dependence of the null vectors can be easily derived by differentiation of 
( p.20| ) and use of the linear spectral equations (p.5|)-( pl6D . This dependence reads: 

\vj) = exp{-A{kj)x - n{kj)t}\vf), (2.21a) 
{v.\ = {vf I exp {A(kj)x + n(kj)t}, (2.21b) 

where \vj^^) and {vf^\ are some constant vectors. 

The vectors in equation (|2.21|) together with the zeros constitute the full set of the 
discrete data necessary to specify the soliton matrix T{k, x, t) and, hence, unique solution to 
the Riemann-Hilbert problem ( p.l| )-( pl2D . Indeed, by constructing the soliton matrix T{k) 
such that the following matrix functions 

0+(A;) = ^+{k)T-^{k), (f)Z\k) = r{k)<^z\k) (2.22) 

are nondegenerate and holomorphic in the domains C+ and C_, respectively, we reduce 
the Riemann-Hilbert problem with zeros to another one without zeros and hence uniquely 
solvable (for details see, for instance, Refs. 0-|^,Q). Below by matrix T{k) we will imply 
the matrix from equation ( p. 22 ) which reduces the Riemann-Hilbert problem (|2.1| )- (|2.2| ) to 



the one without zeros. The corresponding solution to the integrable PDE ( p.7| ) is obtained 
by using the asymptotic expansion of the matrix $(A;) as A; — oo in the linear equation ( p. 51) . 



In the A^-wave interaction model it is given by formula (|2.12|) . The pure soliton solutions 
are obtained by using the rational matrix $ = T{k). 
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III. SOLITON MATRICES FOR MULTIPLE ZEROS 



In this section we consider the sohton solution corresponding to a single multiple zero of 
arbitrary order in the case of an arbitrary matrix dimension N. Such soliton solutions will 
be referred to as the higher-order solitons. We will derive the general formulae for the soliton 
matrices corresponding to an elementary higher-order zero (see the definition below) starting 
from the usual elementary soliton matrices of the Riemann-Hilbert problem. Our formulae 
for the soliton matrices corresponding to an elementary higher-order zero are similar to the 



previously proposed ansatz for the 2x2 Zakharov-Shabat spectral problem |]TT| . However, 
in our approach some essential invariance properties and simple evolution formulae for the 
vector parameters in the soliton matrices are given, which were not known before. Thus we 



simplify the ansatz of Refs. [p!0| , pn^ and put it on the rigorous footing. Although we work 
in the case of involution ( p.4|) , usual for applications in nonlinear physics, our approach is 
valid for the general Riemann-Hilbert problem with zero index. Moreover, we present our 
formulae in a form transferable without any changes to that general case. 

Let and from (|2.1|) each have but one zero of order n, ki and ki, respec- 

tively: 

det = {k- kiY^{k), det ^-_\k) = {k - kiY^{k), (3.1) 

where detip{ki) ^ and det^(fci) ^ 0. The geometric multiplicity of ki (fci) is defined as 
the number of the null vectors in the kernel of $+(fci) ($I^(A;i)), see ( |2.20| ). It can be easily 
shown that the order of a zero is always greater or equal to its geometric multiplicity. It 
is also obvious that the geometric multiplicity of a zero is less than the matrix dimension. 
Before we proceed with the construction of the soliton matrix r(A;) corresponding to the 
multiple zero of order n, two important properties must be pointed out. It is convenient to 
formulate them in the form of two lemmas. 

Lemma 1 Suppose vectors \vj) (1 < j < m) are in the kernel of matrix ^^-(A;!), i.e., 

$+(A;i)|t;,) =0, j = l,...,m, (3.2) 

where m is less or equal to ki's geometric multiplicity. Define the new matrix $+(/c) = 
^{k)^X^^{k) where 

X{k) = I - ^f^P, (3.3) 



m 

P=Y1 \^i)iK-%{vj\, K,, = {v,\v,), (3.4) 

and vectors {vj\ (1 < J < m) are arbitrary but they make matrix K invertible. Then matrix 
^+{k) is also holomorphic in the upper half plane. In addition, if a new vector \w) is in the 
kernel o/$+(fci) and is orthogonal to (TJj| (1 < j < m), i.e.. 
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<^+ih)\w) = 0, {v,\w) = 0, j = l,...,m, (3.5) 

then 

<l>+{h)\w)=0, (3.6) 

i.e., \w) is also in the kernel o/$+(A;i). Furthermore, \w) is linearly independent of\vj) (1 < 
j < m). Similar results exist for matrix $I^(A;) where the multiplication is one the left. 

Remark: it is easy to see that in order for K to be invertible, it is necessary that vectors 
(1 < J < ttl) be hnearly independent. But this condition is not sufficient. However, if 
{vj\ = \vj)'^ (j = 1, . . . , m), then it can be shown that K is invertible. 

Proof. The matrix P is clearly a projector matrix, thus 

Then, expanding the holomorphic function $+(A;) into the Taylor series and recalling equa- 
tion (^.21), we see that 



^^{k)x^\k) = {$+(^0 + {k- k,Y^^ + {k- k,r'^§^ + ...}(!+ ^p) 

= ^+(k,) + ik,-hy^^p+ik,-h)ik-k,y^§^p+..., (3.8) 

which is clearly holomorphic. 
Next, if 

$+(A:i)|w) =0, 

recalling the definition of $+(fc) and expanding $+(A;) into the Taylor series, we get 

<l>+(A;i)|^) + {ki - h /^]^^'^ P\w) = 0. (3.9) 

ak 



Since by assumption, 
thus, 

consequently. 



{vj\w) = 0, j = 1, . . . ,m, 
P\w) = 0, 

^+{ki)\w) = 0. (3.10) 



Lastly, \w) is linearly independent of \vj) (j = 1, . . . , m) because the matrix K is invertible. 
Q.E.D. 

Corollary 1 Suppose the kernel of $+(A;i) is spanned by vectors \vj) {1 < j < m) where 
m is the geometric multiplicity of zero ki. Define matrices x{k), ^+{k) and projector P as 
in Lemma |l|. Then, there exists no vector in the kernel of $+(A;i) which is simultaneously 
orthogonal to {vj\ (1 < j < m). 
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Lemma 2 Suppose that has r independent vectors in the kernel: 



$+(A;i)|i;,) = 0, j = l,...,r, (3.11) 

i.e., rank$+(/i;i) = N — r. Then the following matrix function ^+{k) = ^^{k)x^^{k) , where 
matrix x{k) is as defined in equation l \3.3i ) but with 



P=Y: h){K-%{vj\, K,, = {v,\v,), (3.12) 

has at most r vectors in the kernel at k = ki, i.e., rank$+(A;i) > N — r. Here vectors 
(^jl (1 ^ J ^ ore arbitrary but they make matrix K invertible. 

Proof. This lemma is easy to prove by contradiction. Suppose that there are at least r + 1 
independent vectors \ur+i) in the kernel of $+(A;i) defined above. Then one can 

find a non-zero vector \X) in the kernel of $+(A;i) such that 

(tJ,|X) = 0, j = l,...,r. (3.13) 

Indeed, substitution of the expansion 

\x) = Y.c,\u,) 

into Eq. ( p.l3| ) leads to an under determined, hence, solvable system of equations 

r+1 

^l(Vi\uj)Cj = 0, z = 1, . . . ,r 
i=i 

which have non-zero solutions. But then, according to the second part of Lemma 1, |X) is 
also in the kernel of $+(^1), thus 

\X) = j:C,\v,). (3.14) 

Substituting Eq. ( p.l4|) into ( p.l3|) and recalling that the matrix K is invertible, we find 
that Cj = 0, j = 1, . . . , r, hence X = 0. Thus we have arrived at a contradiction. Q.E.D. 
(Note that a similar lemma is valid for $I^(A;) k = ki with the multiplication on the left.) 

To clarify the implications of Lemma ^ for the soliton matrix T{k) of the higher-order 
zeros, k = ki oi det $+(fc) and k = ki oi det $I^(fc), let us examine the way such matrix is 
constructed. Starting from the solution $±(fc) to the Riemann-Hilbert problem ( p.l| )-( pl2D , 
one looks for the independent null vectors for the matrices $_|_(/ci) and <^Z^{ki): 

$+(A;i)|t;a) = 0, (Ua|<f:'(^i) = 0, z = l,...,Si, (3.15) 

where si is the smaller of ki and fci's geometric multiplicities. Here we allow the two 
geometric multiplicities to be different in general, but they are always the same in the case 
of involution (|2.4|). Next, one constructs the elementary matrix 
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(3.16) 



where 

id 

It can be shown that detxi = (f~f^) ^ ■ If Si < n, where n is the order of the two zeros, then 
one considers the matrix functions ^+{k) = ^^{k)xi^{k) and $I^(A;) = Xi{k)^Z^{k). From 
Lemma |I], we know that matrices $+(/c) and $I^(A;) are also holomorphic in the respective 
half planes of the complex plane. In addition, ki (ki) is still a zero of det$+(A;) (det$I^(A;)). 
Repeating the above steps one gets the elementary matrices Xi{k), Xr{k) such that 
Sl + S2 + . . . + Sr = n. Therefore, 

T{k)=Xr{k)-...-X2{k)xi{k), (3.18) 



where matrices Xi{k) and projectors Pi are as defined in equations (|3.16|) and ( p.l7|) but the 



independent vectors \vii) and {vu] (i = 1, . . . , si) are from the kernels of ($+X]~^ ■. . •■Xil\)(A;i) 
and {xi-i ■ ■ ■ ■ ■ Xi^-^){ki) respectively. 

Lemma |^ indicates that in fact the sequence of ranks of the projectors Pi in the matrix 
T{k) given by equation (|3.18|) , i.e. built in the described way, is non-increasing: 



rankPr ^ rankPr„i < . . . < rank Pi. (3.19) 

This result allows one to classify possible occurrences of a higher-order zero of the Riemann- 
Hilbert problem for arbitrary matrix dimension A^. In general, for zeros of the same order n, 



different sequences of ranks in formula (|3.19|) give different classes of the higher-order soliton 
solutions. In the present paper we consider in detail only the higher-order zeros when the 
sequence of ranks ( p.l9|) is the simplest possible: rankP^ = 1, / = 1, . . . ,n. We introduce 



the following definition. 



Definition 1. In the soliton matrix ( |3. J^ ) corresponding to a higher-order zero ki of a 



Riemann-Hilbert problem, if the ranks of all projectors Pi{l < I < n) are 1, then we call 
this zero an elementary higher-order zero. 

Remark 1: We observe from equation ( |3.19| ) that a higher-order zero (of arbitrary alge- 



braic multiplicity) is elementary if and only if rankPi = 1, i.e., the geometric multiplicity of 
the zero is 1. 

Remark 2: If the matrix dimension N = 2 (as for the nonlinear Schrodinger equation), 
then all higher-order zeros are elementary since rankPi is always equal to 1. 

Below we derive the soliton matrix T[k) and its inverse for an elementary higher-order 
zero. The results are presented in the following lemma. 

Lemma 3 Consider a pair of elementary higher- order zeros of order n: k = ki in C+ and 
k = ki in C_. Then the corresponding soliton matrix T{k) and its inverse can be cast in the 
following form: 
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I + m,...,\Qi))Dik) 



{Pi\ 



\{Pn\ J 



(3.20a) 



^ 1^ 



(3.20b) 



where the matrices D{k) and D{k) are defined as 



( 1 



D{k) 







(fc-fci) 
1 1 



{k-kiY (k-ki) 



V (fc-fci)" ■ ■ ■ (fc-fci)2 (fc-fci) J 



( 1 



D(A:) 



1 



1 \ 



(fc-fci) {k-ki_Y 





(fc-fci)" 



1 1 



(fc-fci) (fe-fcl)2 

n 1 

^ (fc-fci) / 



(3.21) 



and vectors \pj), {pjl \qj) (j 



,n) are independent of k. 



Remark: In ||TT|, the ansatz of the form ( |3.20|) was proposed for higher-order sohtons 
in the nonhnear Schrodinger equation. The lemma above, together with Remark 2 below 
Definition 1, shows that their ansatz is in fact the most general soliton matrix for the 
nonlinear Schrodinger equation. If > 2, their ansatz then is just the soliton matrix for 
elementary higher-order zeros. 

Proof. The representation ( |3.2U| ) can be proved by induction. Consider, for instance, formula 
( p.20a| ). Obviously, this formula is valid for n = 1. In this case, T{k) reduces to an elementary 
matrix x(^)- Now, suppose that formula ( p.20a| ) is valid for n = m. Then we need to show 
that it is valid for n = m + 1 as well. Indeed, denote the soliton matrices for n = m and 
n = m + 1 by r(A;) and T{k) respectively. Then taking into account expression ( |3.18| ) and 
recalling our assumption of the elementary higher-order zero, we have 



T{k)=Xn.+iik)Tik) 



1 + 



k — ki 



m I 



\1j){pl+l-j\ 



(3.22) 



Here, for simplicity of the formulae below, we have normalized the vectors \vm+i) and (Urn+i| 
such that {vm+i\vm+i) = ki — ki. Let us now multiply the two terms in the right-hand side 
of equation ( p. 221 ) and compute the coefficients at the poles: 



k — ki {k — kiY 



+ . . . + 



A 



m+l 



{k-ki)"'+^ 



(3.23) 
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where 



Am+l = \v^+l){Vm+l\A.^ = \Vm+i) {Vm+l\qi) {Pi\, 
A-m = \Vm+l)^m+l\A.m~l + Am 

= \Vm+l){Vm+l\{\q2){Pl\ + I^l)(p2l) + \Qi){Pi\ 
= {\Vni+l){Vm+l\q2) + IQi)) {Pi\ + I Wl) ( 1^1 ) (p2 1 , 

Am-1 = \'Vm+l) {'Vm+l\Am~2 + 

= |^m+l)(Wl|(l^3)(Pll + I^2)(P2| + I^l)(p3l) + l^2)(Pll + {P: 
= {\Vm+i){Vm+l\%) + 1^2)) (Pll + {\Vm+l){Vm+l\q2) + l^l)) {P 
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'2 1 



+ \Vm+l){Vm+l\qi){P3\, 
. . . , 

m 

Al = \Vm+l){Vm+l \ + \qm+l-j){Pj\- 

Define new vectors: 

= \Vni+l){Vm+l\qi), = \Vm+l){Vm+l\qj) + \qj-l), J = 2, 



,m, 



|7„«> = |5„.), (3.24) 

\''^m+l|yi/ 

Then matrices Ai, . . . , A^+i take the following representation: 

I 

Am+2-i = J2\^i+i~j)iPj\^ / = l,...,m + l. (3.25) 



Thus formula ( p.20a| ) is valid for T{k) as well. It is noted that we must also show that the 
denominator in formula ( p. 241) is nonzero. This is easy to show, as {vm+i\qi) is actually (up 
to a factor {ki — fci)™) a product of inner products {vj+i\vj) (1 < j < m), where \vj) and 
{vj\ are the projector vectors in matrix Xj'- 

xAk) = I-^-f^\v,){v,\ (3.26) 
k — ki 



see equations ( p.l8| ) and ( p.22| )]. If {vj^i\vj) = for some j, then Lemma |I] indicates that 



the projector Pj in matrix Xj [see (|3.16| ) to ( |3.18D ] would have rank higher than 1, which 



contradicts our assumption of elementary higher-order zeros. Thus {vj+i\vj) ^ for all j, 
consequently, {vm+i ^ 0. Expression for T'^k) { ^2^ can be proved in the same way. 
Q.E.D. 

In the expressions for r(A;) ( p.20a| ) and r^^(A;) ( |3.20b| ) there are twice as many vectors 



as in the elementary matrices ( p.l6| ) and ( p.l8| ). As the result, only half of the vector 
parameters, namely \pi), . . . , and . . . , are independent. To derive the formulae 
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for the rest of the vector parameters in ( p.20| ) we can use the identity r(A;)r ^{k) = I. The 
poles of r{k)r~^{k) at k = ki, starting from the highest order pole, give: 

r(A;i)bi)(gi| =0, 

1 dr(A;i 



m){\P2){qi\ + \pi){q2\) + j^^-^\pi){qi\ 



m) (b3)(gi| + \P2){q2\ + bi)(g3|) + ^^^^ {\P2){qi\ + \Pi){q2\) + = 0' 



Hence, we obtain: 



r(fci)bi) = o, 
1 

1! 

1 dr(A: 



m)b2) + ^^bi) = 0. 



1 d2r(fci) 



r(fci)b3) + ^^^^h) + ij^^m^lPi) = 0, 



1! dk 



2! dF 



r(A;i)b„) + 



1 dr(A; 



1! dk 



+ ...+ 



1 d'^-T(A; 



>-!)! dk 



(3.27a) 
(3.27b) 

(3.27c) 
(3.27d) 



Equations ( |3.27|) can be written in a compact form for the following matrix T{k): 



r(fci 



^ IP.) ^ 



y \Pn) J 



( 



V{k) 



J. Ar 

1! dfc 



1 d"""" -p F 
\ (n-1)! dfe"-!-"- ■■• lldfe-"- / 





r 



1 d 



(3.28) 



Note that, as a block matrix, r(A;) has (lower-triangular) Toeplitz form, i.e. along each 
diagonal it has the same (matrix) element. 

In much the same way, by considering the poles at = fci in Y{k)Y'^{k\ one derives the 
following formula 



((Pil,---,{P„l)r(*:i) = 0, T(k) 



1 Af-1 

1! dfc 

r-1 



1 d" 



(n-1)! dfc" 



1 Af-1 

1! dfc-^ 



V 











. (3.29) 



Equations (|3.28|) and ( p.29|) allow us to find the expressions for the dependent vector 
parameters. For convenience of the presentation, let us introduce the following fc-dependent 
vectors: 
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Then, by reordering the summation in (|3.20|) we get 



r(A;) = J + (|gJ,...,|gi)) 



T~\k) = I+{\Z,,{k)),...^Z,{k))) 



Let us now substitute the expression ( |3.31a) ) into equation ( p. 281 ) and solve for \qxi-, 
We have 



(^i(fci)bi) 



b2) + (|gn),...,|gi)) 



(Zi(A;0b2) + M(^i(fci)bi) 



b3) + (|g„),...,|gi)) 



V 



Hence 



where 



/ n_ \ 



(Zi(fci)bi) (^i(A:i)b2) + ^ii(^i(fci)bi) ... g („io!dV-' (^i(^i)bO 
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Similarly, we get 



where 



' in ^ 



\ {Pn\ J 



(3.34) 



{Pl\Zn{kl)) 



{Pi\Zi{ki)) 



(P2|^n(A:i)) + ^ii(Pi|Z„(A;i)) ... (p2|Zi(A:i)) + i,^(Pi|Zi(fci)) 



v 



E 



{n-l)\ dk"-l {Pl\Zn{kl)) 



E 

1=1 



1 d" 



{n-l)\ dfe'^ 



'-{pi\Zi{ki)) 



J 



(3.35) 



In terms of the independent vector parameters, the soliton matrices (|3.20a ) and ( 3.20b ) 
can be rewritten as 



T{k)=I-{\p,),...,\pn))lC-^D{k) 



(3.36) 



T-\k) 



(bi),...,|p„))D(A:)/C-i 



I 



(3.37) 



where matrices /C and /C are given in equations ( p.35| ) and (|3.33|) . 

The soliton matrices given by ( |3.36| ) and ( p.37| ) possess invariance properties. The in- 
variance is the transformation of the independent vector parameters which preserves the 
form of the soliton matrices and equations defining the vector parameters, i.e. equations 
( p.28| )-( pl2^ ). Let us first consider transformations of vectors \pj) (j = 1, . . . ,n). Suppose 
these vectors are transformed as 



(bi),---,bn)) = (|pi),...,|p„))5. 



(3.38) 



where i? is a /c-independent matrix which, in general, depends on (x,t). Here the vectors 
(Pjl (j = 1) • • • ! remain intact. Simple calculations show that the new vectors . . . , 
satisfy equation ( p.28|) if and only if the matrix B has upper-triangular Toeplitz form. 
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B 



b2 
hi 




(3.39) 



V 



h2 
hi 



I 



Further, we note that under the transformation ( |3.38| )-( p.39| ) the matrix /C transforms as 



(3.40) 



where matrix /C is as given by equation (|3.33|) but with vectors \pj) replaced by the new 
vectors \'Pj)- From formulae ( 3.38|) and ( p.40 ) it is seen that the form ( p.36 ) of matrix r(A;) 
is preserved. We still need to show that for matrix B of the form (|3.39|) , the transformation 
( p.38| ) also preserves the form ( |3.37| ) of matrix r^^(A;). Notice that matrix -D(fc) also has 
upper-triangular Toeplitz form, thus -D(fc) and B are commutable. Utilizing this property, 
we can easily show that under the transformation ( |3.38|) , matrix /C transforms as 



(3.41) 



where /C is given by equation (|3.35|) but with \pj) replaced by \pj) ■ Thus the form of matrix 
r~^(/c) is also preserved. In short, soliton matrices ( p.36| ) and (|3.37| ) are invariant under the 
transformation ( |3.38| ) with matrix B given by ( p.39| ). 



Similarly, we can show that soliton matrices ( p.36|) and ( |3.37| ) are also invariant under 
the transformation 



\Vi\ 



B 



\Vi\ 



\ kPn\ ) 



(3.42) 



where the /c-independent matrix B (which, in general, depends on (x, t)) has lower-triangular 
Toeplitz form. 



B 



^ hi 



62 hi 







y hn 



hi 



(3.43) 



and vectors \pj){j = 1, . . . , n) remain intact. 
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Summarizing, we conclude that soliton matrices ( p.36|) and ( |3.37|) are invariant under 
the triangular Toeplitz transformations 



{Pl\ 



\{Pn\ j 



B 



{Pi\ 



\ {Pn\ J 



(3.44) 



of the independent vectors \pj) and {pj\ (1 < j < n). Here B and B are arbitrary lower and 
upper triangular Toeplitz matrices, respectively, in general (x, t)-dependent. Here we point 
out that the invariance transformation found in [|ll| is given by bj = bj = {2 < j < n — 1) , 



6, 

i.e., only and 6„ being non-zero. Thus it is just a special case of the invariance 

property of the soliton matrices. 

The invariance transformations indicate that arbitrary sets of vectors \pi), ■ ■ ■ , \pn) and 
{pi\, . . . , satisfying equations (|3.28| ) and (|3.29| ) can be chosen as the independent vector 
parameters. This is, in fact, also a necessary condition for such vector parameterization of 
the soliton matrix to be self-consistent. 

Now let us derive the (x, i:)-dependence of the vector parameters which enter the soliton 
matrix. We can start with the fact that the soliton matrix r{k,x,t) must satisfy equations 
^SD-fP.GI) with some potentials U{k,x,t) and V{k,x,t): 



d^T{k, X, t) = T{k, X, t)A{k) + U{k, x, t)T{k, x, t), 
dtT{k, X, t) = r{k, X, t)Q{k) + V{k, x, t)r{k, x, t). 



(3.45a) 
(3.45b) 



The derivation is based on the use of equations ( p.28| ) and ( p.29| ) (quite similar to the 
derivation of equations ( |2.21| ) in section y). First of all we need to find the equations for the 
triangular block Toeplitz matrices T and T. To this goal one needs to differentiate equations 
( ^.45| ) with respect to k up to the {n — l)-th order. It is easy to see that, for instance, the 
equations for the F have the same form as equations ( |3.45| ): 



d^r{k, X, t) = T{k, X, t)A{k) + U(A;, x, t)T{k, x, t), 
dtVik, X, t) = r{k, X, t)n{k) + V(A;, x, t)T{k, x, t), 

if we introduce the lower-triangular block Toeplitz matrices A, O, U, and V: 



(3.46a) 
(3.46b) 



/ 



A = 



A 



1! dfc^^ 



A 

1 d 



\ (n-1)! dfc 



"-1^^ • • • 1! dfc^^ 




A 



n 



J__d 

1! dk 



1 d 



1 d"~" n -L-S-O 

\ (n-l)! dfc"-i" • • • 1! dfe" 







(3.47) 
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u 



u 

J__d_ 

1! dfc 



U 



\ (n-iy. dfc" 



u 

1! dfc^ 





u 



V 



V 

1! dfc 



V 



1 d"-i T/ 
\ (n-1)! dfc"-i 



J__d_ 

1! dfc 



V 




V 



(3.48) 



Indeed, this is due to the fact that the matrix multiphcation in ( |3.46| ) exactly reproduces 
the Leibniz rule for higher-order derivatives of a product. Similarly, using the equations for 
r~^, one finds that 



d^T{k, X, t) = -A{k)T{k, x, t) - T{k, x, t)V{k, x, t), 
dtT{k,x,t) = -U{k)T{k,x,t) -T{k,x,t)Y{k,x,t), 

for the upper-triangular block Toeplitz matrices A, f2, U, and V: 



A 



1! dfc^^ • • • (n-1)! dfc 



—A 



A 



1! dfc^*- 

A 



n 



" 1! dfc" • ■ ■ (n-1)! dfc 



n-1 ^ ' 



n 



1! dfc" 



(3.49a) 
(3.49b) 



(3.50) 



U 



u 



1 d"-l_fj \ 



(n-1)! dfc 











1! dfc^ 

u 



J 



V 



V 



1 d"-l_y \ 



(n-1)! dfc 



1! dfc ^ 

V 



... 



(3.51) 



The {x, t)-dependence of the vector parameters \pi), . . . , |p„) and {p-^\, . . . , {p^\ can be found 
by differentiation of equations ( |3.28| ) and ( p.29| ) with the help of equations ( |3.46| ) and ( p. 491 ). 
First, we note that for commuting matrices the corresponding block Toeplitz matrices as 
introduced above also commute with each other. Second, it is shown in the Appendix 
that for a diagonal matrix [e.g. A{k)x + Q{k)t] the operation of raising to the exponent 
commutes with the construction of the block Toeplitz matrix. Therefore, taking into account 
the invariance property, we find the (x, t)-dependence of the vector parameters as 



bi) 



y \Pn) J 



exp {— A(fci)x — ft{ki)t} 



(3.52a) 
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((Pil, • • • , (Pnl) = • • • , {P^^D exp {A{h)x + n{h)t} . 



(3.52b) 



Here the superscript "0" is used to denote constant vectors and the exponents stand for the 
triangular block Toeplitz matrices: 



/ 



exp {—A{ki)x — Q{ki)t} 



E{h) 



V (n-l) 











E{ki) 



(3.53a) 



exp 



E~\k,) ^,f,E-\k,) 
E-^{ki) 



\ 



1 d"-i T?-i(i: \ \ 

(n-l)! dfc"-!-^ '^"'1^ 



h.i!:E-\k,] 
E'\k^) 



I 



(3.53b) 



where E{}z) = exp {—A{k)x — Q{k)t}. After the temporal and spatial evolutions for vectors 
\pj) and {pj\ have been obtained as above, the corresponding higher-order soliton solution 
can be obtained from equations (^), ( ^TT^ ), (gM^) and {^^. 



IV. APPLICATION TO THE THREE WAVE INTERACTION MODEL 



Here we apply the theory developed in the previous section to the three-wave interaction 
model (|2.14|) . The three- wave model has wide applications in nonlinear physics. For instance, 
under the additional constraint Uj = iqj where qj are real variables, it describes the "exact 
resonance" in parametric interaction of three wave packets, while under the reduction of the 
dispersion laws ( |2.11| ): = —ai, a2 = 0, 63 = —61, 62 = and the condition U2 = —Ui, it 
models the generation of second harmonics. The usual (fundamental) soliton solutions to 
the three- wave interaction model have been well studied (consult Ref. 0). Such solitons 
approach sech profiles as t ±00 on the characteristics x — Vjt = const. 

Let us consider the simplest higher-order solitons in the three-wave system: solitons 
which correspond to an elementary higher-order zero of order 2. Here we should take into 
account the involution property given by equation (|2.4|). For instance, we have 



k = k* 



{Pj\ = IPj)^ 



(here and below the overline is associated with the Hermitian or, in the case of scalar 
quantities, complex conjugation). Then the soliton matrix reads 
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{\pi)Ap2))k: 



{Z,{k)\ 



(4.1) 



where 



V 



{Z,ih)\p,) {Z,ih)\p2) + i{z,ih)\p,) 



(4.2) 



and 



[Z2m 



{pi\ 



{P2 



{Pl\ 



k — ki ^ " k — ki {k — kiY 
The (x, t)-dependence of the vector parameters \pi), \p2) has the following form 

IpP) 



\P1) 

\P2) 



\ \P2l J 



E{ki) 
y ^^E{kr) E{k^] 



\ ( 1^(0)) \ 



V 



E{k, 



-iki(Ax+Bt) 



(4.3) 



We denote fci = ^ + ir], where ^ and rj are real numbers (r^ > since ki hes the upper 
half plane of the complex plane), and choose the following parameterization of the constant 
vectors \pf^) and \p^2^): 



\pf^) = 2ir] 



( 0(1) \ 



/ 0(2) \ 

n(2) 
V^3 J 



(4.4) 



where ^j*'*'s are complex constants. It is noted that due to the invariance property ( ^.44p , 
where the matrix B contains two arbitrary complex constants, we have 2 free components 
in each vector in formula ( [4.4|) . Hence there are 10 free real parameters (including ^ and rj) 
in the higher-order soliton solution. 

The (x, t)-dependence of the components of the vector parameters reads 



Pij = 2ir]9j e 



(1) f,/2-ix, 



P2j 



Jj/2-iXj 



where 



fj = 2r]{ajX + hjt), Xj = ^{^jX + bjt), j = 1, 2, 3. 
By simple calculations we obtain the elements of matrix K. as 

= -2^r^E l^f |V^ ICu = -j:(0f^Of + (/, - 1)|0«I^^ 



(4.5) 



(4.6) 



(4.7) 
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j=i ^ ^ i=i 

(4.8) 

It is easy to verify that the determinant of /C is 

detZ=-± (l^.W^f V + he^ef^ - (^ef + (/. - fM''0f^') e^^-^ (4.9) 

which is always non-zero: 

For the sohton solution corresponding to the matrix (f4. 1|) we need the first-order term 
of its asymptotics as k ^ oo: 

r(^) = (^22bl)(Pl| +^llb2)(p2l -^12bl)(p2l -^2lb2)(Pl|) . (4.10) 

Using formulae ( [4.5|) for \pi) and \p2) and ( |4.7| , ^78|) for the elements of /C we get 



det:^ ^ti ' ' ^ ^ 



where 



Imj 



efe^"^ (/, - - 1) + - ef^e^p\ IC^}" - /. + 1) + o^W - 



-^W^Vf^P- (4.12) 



The three nonlinear waves Ui,U2,U3 are given by formula (|2.13|) . Thus 



ui = Voi -02 r^a^ M2 = Va2 - OsT^a^ % = V«i - «3 T^g^ (4.13) 



where r^-j'' are given by Eqs. (^4.11|) and ([4.12|) . To be explicit, our soliton solution corre- 



sponding to an elementary higher-order zero of order 2 in the three-wave interaction system 
is 



^!!IV^lZ^^^h+h)/2-^ix^-X2) ^ C'l2,e^^ (4.14) 



Ml 

det /C 



^!W«^a^g(/,+/3)/2-.fe-x3) ^ C23,ef\ (4.15) 



^2 

det IC -^^ 



'^^!lV^lZ^^ih+h)/2-^(xr-xs) J2 Cl3,e^^ (4.16) 



^ det /C j-^^ 
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where det /C, Cijk, fk and Xk are given by Eqs. (^^) , (|4.9| ) and (|4.12|) 



The above solutions are fairly complicated. But some information about them can be 
gained from considering the asymptotics as t — > ±00. Evidently, the t-asymptotics is nonzero 
only on the characteristics: 

zi = (/i - /2)/2 = ?7(ai - a2){x - Vit), 
Z2 = (/2 - /3)/2 = ri{a2 - az){x - V2t), 



z-i = (/i - /3)/2 = r]{ai - a3)(x - V3t). 

The asymptotic formulae depend on the relation between the velocities of the waves. For 
definiteness, let us choose V2 < fi. This is equivalent to the condition e > in view that 



- {V1-V2), 



ai — 03 



as follows from formulae (|2.15|) - (|2.16|) . Then the condition (|2.17|) requires that lies between 
V2 and fi: 

V2 <V3< Vi. (4.17) 



Further, we notice that any solution ui,U2,U3 of the three-wave interaction model (p.l4|) 



in the case of the opposite inequality, i.e. fi < ^3 < V2, is mapped onto the solu- 
tion satisfying the inequality ( [4. IT] ) by the following transformation: Uj{x,t]Vi,V2,V3) = 
—Uj{x,—t;—vi,—V2,—V3). Thus, the case of vi < V2 is easy to recover (it describes the 
reverse process to that of V2 < fi). 

The asymptotic formulae also depend on whether some of the components in vectors O^^"* 
and 6^'^^ are zero or not. We first consider the generic case when none of the parameters ^j^^ 
for j = 1, 2, 3 is zero. Define the following real quantities: 



1 I I 7/1 I I " it TYh I 



1 /^^p 


0(2) 
m 


2 [el^^ 








V 



and denote 



V^S^^ = arg(#), ^,^ = ^\^>-^(^^---:^, (4.18) 



^12 = ^1- "12, ^23 = ^2 -"23, ^l3 = 2;3-ai3- (4.19) 



Then, simple calculations show that the asymptotics of the waves ( [4.14| ), ([4.15|) and (f4. 16|) 
are as follows: 

Ui ^ 0, U2 ^ 0, t ^ — oo; M3 — > 0, t ^ 00, (4.20) 



23 



U2 



■ 0-2 



(^12 + gi2) sinh zx2 - (1 + i(y\2) cosh Zx2 Mv^2-iiz^2l'n) 



cosh^ Zi2 + (2:12 + ^'12) 



C^12 



t — > 00, 



(4.21) 



Ms 



cosh Z23 + (^23 + ^'23)^ + 0-23 

^ (^13 + gia) sinh - (1 + ^^13) COshZis ^(^,3-£.l3/ry) 

^"/V "1 "3 1 2 , / , \9 , 2 ' 

cosh 2i3 + (Zi3 + ^13)^ + 0-13 



t 00, (4.22) 



t -00. (4.23) 



We see that as t — —00, only the pumping wave U3 is non-zero, while as t ^ cxo only 
the elementary waves ui and U2 are nonzero. Thus in the generic case of higher-order 
solitons under condition ( |4.17| ), the solution describes the breakdown of the pumping higher- 
order soliton M3 into the higher-order solitons of the elementary waves Ui and U2- [For the 
opposite inequalities in formula ( [4. IT] ), the solution describes the reverse process: merger 
of the two elementary waves ui and U2 into the pumping wave M3.] These properties are 
identical to fundamental solitons (see, for instance, pp. 174-184 in Ref. 0). However, 
differences between higher-order solitons and fundamental solitons are also obvious: none 
of the asymptotics ( [4.21| ) to ( [4.23| ) of the higher-order solitons is sech-shaped, while the 
asymptotics of fundamental sohtons are all sech-shaped. 

Asymptotics (|4.20| ) to ( [4.23| ) are invalid in the non-generic cases when at least one of the 
parameters 6j^\ j = 1, 2, 3, is zero. Consider first the case of ^i^^ = and 6^2^ and 6^3^ bemg 
non-zero. There are two possibilities depending on whether 61 is zero or not. 

(a) If 9^^ 7^ 0, then the asymptotics of the waves ([4.14|) , ([4.15|) and ( [4.16| ) become 



Ul 



0, 



t — i> —00: 



U3 



0, t 



00; 



(4.24) 



Ul -ir]^ai - a2e*^^^''"^^''"^"^/^^sech(2i - /3i), t 00, 



(4.25) 



U2 ^ -ir]^a2 - a3e*('^2''-4''-«-2/r,)ggpi^(^^ _ ^^^^ ^ 



-00, 



(4.26) 



U2 22771/02 - ^3 



(^23 + g23) sinh Z23 - (1 + ^g"23) COsh Z23 ^ 

cosh^ 2:23 + (2:23 + Q2zY + <yh 



,ii'P23-(.Z23/v) 



t 00, (4.27) 



Us -ivVai - a3e'^'^^''"'^^''"^"^/''^sech(z3 - /S^), t 
where parameters Pj (j = 1, 2, 3) are defined as 




/?2 




/?3 



-00, 




(4.28) 



(4.29) 
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The above asymptotics have two important features. One is that as t — > — cx), both and 
U2 waves are non-zero. Thus, the higher-order sohton solution ( [4.14| ) to ([4.16|) in this non- 
generic case does not describe the breakdown of the pumping wave u^. Instead, it describes 
a new process: 



U2 + Us Ui + U2. 



(4.30) 



This is very different from fundamental solitons. The other feature is that the waves U2 and 
M3 as t ^ — 00 and the wave ui as t ^ 00 all have sech profiles, but the wave U2 as t ^ 00 
is the higher-order soliton. 

(b) li Oi = 0, then the higher-order solitons ( |4.14|) to ([4.1 6| ) are degenerate: 



(4.31) 



Ml(x, t) = U3{x, t) = U2{x, t) = U2'\x — ^2^), 

where U2'\x) is the initial solution of U2- This is a trivial solution. 

If two components of the vector 9^^^ are zero, then higher-order soliton solutions (|4.14|) 
to ( [4.16[ ) reduce to fundamental-soliton solutions or trivial solutions. For instance, if 6i = 



62^ = 0, 6*3^'' 7^ 0, 9i'^ 7^ 0, and 62'^ 7^ 0, then the asymptotics of the waves become 



Ui 



0, U2 0, t ^ —00; 



M3 



0, t 



00; 



(4.32) 



ir/Vai - a2e'('^i"-'^2"-?-i/^)sech(;zi - /3i), t 00, 



Ml 



(4.33) 



U2 ^ -ir]^a2 - a3e*(^2''~^i''-C-2/^)sech(z2 - h). t 00, 



(4.34) 



Ms 



, ., (2) (1) , ; N 



sech(2;3 — [3^), t —>■ —00, 



where f]^ is defined in equation ( [4.29| ) and 
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This is the fundamental soliton solution. If one or both of ^ and ^^'^ is zero, the solution 



3(2) 



(4.35) 



(4.36) 



is trivial (similar to (|4.31|) ). We note that when 6^^ = 6^^ = 0, then ^3^^ can not be zero. 



because otherwise, the denominator det /C in the solution is zero. 

It turns out that consideration of the case when 6^^^ = is similar to the above case of 
6^^ = with the only difference that now the elementary waves Ui and U2 are interchanged. 



For instance, when 6^^ = 0, OY^' 7^ 0, 62'^' 7^ 0, and ^3^'' 7^ we have the following (also new) 
process: 

U1 + U3 ui + U2, (4.37) 
where the waves Ui and M3 as t ^ —00 and the wave M2 as t — > 00 all have sech profiles, 



3(2) 



while the wave Ui as t — cxd is the higher-order soliton [see equation (|4.22| )]. 
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The only (different) case which is left to consider is the case of = with e^^^ ^ 
and Ol^^ 7^ 0. The asymptotics depends on whether 6^^ is zero or not. In the former case, 
i.e. 62 = 0, we have a degenerate solution, 

ui{x, t) = U2{x, t) = 0, us{x, t) = uf\x — Vst), 



which is similar to solution ([4.31|) . If however 62'^ 7^ 0, then the asymptotics of the waves 
(|;T|)-(|;Tg) are as follows: 



Ml ^ 0, U2^0, -00, (4.38) 

.3 ^ 2^,V^^s ^"'' ^ ''i 7^.13 ^,,.3-,.W.), t^_^. (4.39) 

cosh zis + {zi3 + gi3y + crf3 

Ml ^ -ir]Vai - a2e*('^i'^-'^2''-«^i/^)sech(2i - /3i), t ^ 00, (4.40) 

U2 -tvVa2 - a3e*(^2''-4''-e^2/^)sech(z2 - t 00, (4.41) 

U3 -ivVai - a3e*^^i'^"^3''-5^3/r;)gg^j^(^^ _ ^ (4_42) 

where is defined in equation ([4.36|) and 

These asymptotic formulae describe yet another new process: 

^^3 Ml + M2 + Us, (4.44) 

where waves ui, U2, and M3 as t — > 00 all have sech profiles, while the pumping wave as 
t — > — 00 is more complicated. Thus, this process describes a breakup of the higher-order 
pumping wave into three sech waves, the two elementary waves and the pumping wave. 

Lastly, we present the graphical pictures of the above higher-order solitons for both 
the generic and non-generic cases. In all figures, the common solution parameters are 
(ai, 02, as) = (2, 1, -1), (61, 62, h) = (-0.5, 2, 1), ^ = 1, = 1, and ^(2) = 1 + ^, 2). 
Only the vector 9^^^ is different. It is easy to check that for these parameters, the inequal- 
ity (|4.17|) holds, thus the asymptotics of these higher-order solitons have been described in 



^^ = in^h f^^ = ^^\\±yJ- (4-43) 



the previous text. In all figures, the solid lines are |ui|, the dashed lines are 1^21, and the 
dashed-dotted lines are \u3\. 

First, we illustrate the generic solution with 6^^^ = (1,-i,— 1) in Fig. 1. As we can see 
from this figure as well as the asymptotics ( [4.20| ) to ( |4.23| ), as t —00, only the pumping 



^3 solution is non-zero. As t ^ 00, this M3 wave breaks into elementary ui and U2 waves. 
This process is similar to fundamental solitons. But there is a difference: the asymptotics 



26 



of each wave in Fig. 1 has a complex structure which signals that it is a higher-order 
soliton instead of a fundamental soliton. Next, we let 9i^^ approach zero. Specifically, we 
let e^^^ = 10-^ while the o!^^ and 6^^^ values remain the same. The corresponding soliton 
solution is illustrated in Fig. 2. We see that in this case, the pumping wave at t — > — oo 
splits into two sech pulses. As time moves on, the front sech pulse breaks into Ui and 
U2 sech pulses. Then this U2 sech pulse and the back sech pulse interact. The final 
outcome is two ui sech pulses moving in the positive x direction, and a higher-order U2 
wave moving in the negative x direction. Thirdly, we consider the non-generic case where 
6i^^ = 0, while 62^^ and ^3^'' still do not change. This soliton solution is illustrated in Fig. 
3. We see that as t ^ —00, both the M3 and U2 waves are non-zero and sech-shaped. After 
their interaction, the pumping M3 wave is depleted, and a new ui sech wave and a higher- 
order U2 wave are created. We note that this M2 + M3 — > Ui + U2 process is novel, and it 
has not been carefully investigated before. Fourthly, we consider the non-generic case where 
6[^^ = 9^2^ = 0, and 6*3^'* is still —1. This solution is illustrated in Fig. 4. We see that it is the 
same as a fundamental soliton solution, and it describes the process of a pumping U3 sech 
wave breaking into two elementary ui and U2 sech waves. Thus, our higher-order soliton 
solution reduces to a fundamental soliton solution as a special case. Lastly, we consider the 
non-generic case where 62^^ = while 6^^ = 1 and ^3^^ = —1 as in Fig. 1. This solution is 
shown in Fig. 5. As we can see, as t ^ —00, the only non-zero wave is the pumping wave 
M3, which is a higher-order soliton. As t — * 00, this pumping wave breaks up into a sech 
waves in each component. Thus, this is the new ^ ui + U2 + U3 process which we have 
presented in the text above. 

We conclude this section with some comments on the soliton solutions to the three-wave 
model corresponding to the higher-order zeros of order n > 2. If the higher-order zero is 
elementary, i.e. when the sequence of ranks in formula (|3.19|) is rankP,- = 1, j = l,...,n, 
the corresponding soliton solutions can be derived using the soliton matrix ( |3.36| ). However, 
there are two other possible sequences of ranks in formula ( 3.19| ), namely 



(a) rankPj- = 2, j = 1, ...,r, n = 2r; 

(b) rankPj = 2, j = 1, ...,r; rankP,- = 1, j = r + 1, ...,r + s, n = 2r + s. 

We note that the soliton matrix for the sequence of ranks (a), which corresponds to 
the higher-order zero of order 2r, has an equivalent soliton matrix corresponding to the 
elementary higher-order zero of order r. Indeed, let us consider the soliton matrix in the 
representation ( |3.18| ), where each Xj(^) is defined similar to formula ( |3.16D [Pj substituted 
for Pi] with rankPj- = 2. Consider the following procedure. First, define new projectors 
Qj = I — Pj. Evidently rankQj = 1. Second, multiply the soliton matrix T{k) (|3.18|) by a 
scalar quotient, 

such that each Xji^)y j = 1) ■■■y^y gets a multiplier {k — ki)/{k — ki). We have 
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and 



f{k) = Xr{k)...Uk)Uk). (4.46) 



Evidently, the new matrix r(A;) in (|4.46| ) satisfies the hnear system of equations (|3.45a|) - 



( P .45131) for the original matrix r(A;). Furthermore, it corresponds to an elementary higher- 
order zero of order r, though now in the complementary half plane: k = ki. It is noted that 
in this section we considered a soliton matrix r(A;) corresponding to a zero ki = ^ + 17] lying 
in the upper half plane, i.e. with 77 > 0. However, the case of 77 < is admissible as well. The 
only significant change would be in the asymptotic formulae, and the effect of this change 
is similar to reversing the time variable: t —t. Thus, the sequence of ranks in case (a) 
brings no new higher-order soliton solutions as compared to the simple sequence of ranks, 
but the solution process is reversed. For the fundamental soliton solutions, a similar fact has 
been noted in Ref. 0, where it is mentioned that the fundamental soliton corresponding to 
the projector of rank 2 describes the three-wave interaction process which is reverse to that 
of the soliton solution corresponding to the projector of rank 1. 

There is no transformation similar to ( [4.45| ) for case (b) (similar multiplication will 



produce a rational matrix function having poles in both half planes, thus such a matrix 
does not belong to the class of soliton matrices). Higher-order soliton solutions in this case 
require construction of the soliton matrices for non-elementary higher-order zeros and will 
be addressed in a forthcoming paper. 



V. CONCLUSION 

We have proposed a unified and systematic approach to study the higher-order soliton 
solutions of nonlinear PDEs integrable by the Riemann-Hilbert problem of arbitrary matrix 
dimension. We have derived the soliton dressing matrix for the elementary higher-order zeros 
in the N x A^-dimensional spectral problem, i.e., zeros having the geometric multiplicity 1. 
The associated higher-order solitons in the A^-wave system have also been obtained. We have 



also clarified that the soliton dressing ansatz proposed in ITT] is the general soliton matrix 



for the nonlinear Schrodinger equation (where N = 2), thus the soliton solutions obtained 



in [11] are the most general higher-order solitons in the the nonlinear Schrodinger equation. 
For N X A^-dimensional spectral problems the soliton dressing ansatz of |11| corresponds to 
the elementary higher-order solitons. 

We have applied our theory to the three-wave interaction model, and the simplest higher- 
order soliton solution has been obtained. The generic case of this solution describes the 
process M3 ^ ui + U2, similar to fundamental solitons. But each wave involved here is 
higher-order. The non-generic case of this solution could describe three new processes. The 
first two are similar to each other: 
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Ml + Ms ^ Ml + M2, 
M2 + Ms ^ Ml + U2. 



Here the waves on the left are all sech waves; the waves on the right are a sech wave and a 
higher-order wave. The third process reads 

Ms ^ Ml + ^2 + Ms, 

where the pumping wave on the left is a higher-order wave, and the waves on the right all 
have sech shape. The non-generic solutions could also reduce to fundamental solitons or 
trivial solutions as special cases. 

We anticipate that the higher-order soliton solutions will have wide applications. First 
of all, the new processes they describe may find physical applications where three-wave 
interaction takes place. Second, as it has been mentioned in Ref. |T0|, the higher-order 



soliton solution describes a weak bound state of solitons, thus it may appear in the study 
of the train propagation of solitons with nearly equal amplitudes and velocities in nonlinear 
integrable PDEs. The usual approach in the analytical study of the soliton trains is reduction 
of the governing equations for the soliton parameters to the complex Toda chain (consult, 
for instance, Refs. ||28|-pl|). The higher-order soliton approach may provide an alternative 
to this study. Thirdly, multi-hump solitary waves in the non-integrable nonlinear PDEs 
can be another field of application of the higher-order solitons. For instance, the so-called 
multisoliton complexes, or more precisely, oscillatory and stationary solitons observed in 
an oscillating water trough [P^-pB| and subsequently reproduced in numerical simulations 
3^-p6[ of the governing parametrically driven, damped NLS equation may have the same 



relation to the higher-order solitons as the usual solitary-wave solutions of the non-integrable 
PDEs to the fundamental solitons. Analytical study of the soliton complexes needs the 
perturbation theory for the higher-order solitons, just as the study of usual solitary- wave 
solutions needs the perturbation theory for the fundamental solitons. The perturbation 
theory for the higher-order solitons can be developed in a similar way as it is done for the 
fundamental solitons (see for instance Ref. p7|-|39|]). Such a theory is left for future studies. 

Lastly, we point out that the soliton matrices for the elementary zeros serve as the 
building blocks for the general case of zeros with arbitrary geometric multiplicity. This 
work is in progress and will be reported in a forthcoming paper. There the most general 
higher-order soliton solution for the A^-wave system will be given. 
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APPENDIX A: EXPONENT OF THE TOEPLITZ MATRICES 



Here we show that for a diagonal matrix M{k) the exponent of the block Toeplitz matrix, 
defined as in formula ( |3.47| ), and the block Toeplitz matrix of the exponent of M{k) coincide. 
As the derivatives -^^M{k) commute with each other, it is enough to prove this statement 
for a scalar function. Consider, for example, the lower-triangular Toeplitz matrix of a scalar 
function f{k): 



I 



f 











_d_ 
dk 



f f 



n — lf AUn-2 f • • • f 



(Al) 



1 ' Lm 



(A2) 



\ dk" 

It can be rewritten in the following form 

df d"-V 

Note the product rule for the "diagonals" : HjHi = Hj+i and that for j + i > n — 1 the 
product is zero. Therefore, the exponent of F is a finite sum of the diagonals Hj: 

exp(F) = coHo + CiHi + ... + c„_ii/„_i, (A3) 

where Cq, . . . , c„_i are constants. Due to the formula Hj = Hi, j = 0, . . . ,n — 1, computing 
the coefficients Cj is equivalent to taking the finite sum of the first n terms of the Taylor 
expansion of an equivalent scalar function: 



exp < 2^ - - . 



1=0 



dk^ 



Co + cie + . . . + c„_ie"-^ + C{e"} 



(A4) 



where e is the parameter of the Taylor expansion which represents Hi. On the other hand, 
computing the Taylor expansion reduces to taking derivatives with respect to k of exp{/(fc + 
e)} at e = 0: 



,=0 



dk^ 



exp{/(A; + e)} + 0{e"} 



exp{/(.)} + ^^exp{/(.)}e + . . . + ^^.^ ^ u ... .n-. 



-exp{/(A;)}e"-i + 0{en. (A5) 



Therefore 



thus. 



1 d^' 
j! dk^ 



exp{/(/i;)}, j = l,...,n-l. 



(A6) 



exp{F} = ifoexp{/(fc)} + if^l A exp{/(A;)} + . . . + H^^i ^ f^^\ e^^{f{k)}. (A7) 



Q.E.D. 
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FIG. 1. A generic higher-order sohton solution which describes the breaking of the higher-order 
pumping Us wave into higher-order elementary ui and U2 waves, i.e., the M3 — > «i + U2 process. 
Here, the solution parameters are (01,02,03) = (2,1,-1), (61,62,63) = (—0.5,2,1), ^ = 1,7/ = 1, 
0(1) = (1, ^, -1) and ^(2) = 1 + 2). In all fi gures here and below, solid lines are |ui|, dashed 
lines are |tt2|, and dash-dotted lines are |u3|. 
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FIG. 2. Another generic higher-order soliton solution with a very small value. Here 



01^^ = 10 ^, while the other solution parameters are the same as in Fig. 1. 
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FIG. 3. A non-generic higher-order sohton solution which describes the U2 + U3 ^ ui + U2 



process. The solution parameters are the same as in Fig. 1 except that 6^^ = now. 




FIG. 4. Another non-generic higher-order soliton solution which describes the breaking of the 
Us sech wave into ui and U2 sech waves. The solution parameters are the same as in Fig. 1 except 
that e[^^ = 9^2^ = here. 
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FIG. 5. A non-generic higher-order soliton solution which describes the U3 ^ ui + U2 + U3 
process. The solution parameters are the same as in Fig. 1 except that 6^^ = here. 



34 



